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^ . Abstract 

We prove a Mahoux-Mehta-type theorem for finite-volume partition functions of SU{Nc > 3) 
gauge theories coupled to fermions in the fundamental representation. The large-volume limit is 
^\f ' taken with the constraint V^ << 1/m^. The theorem allows one to express any k-point correlation 

function of the microscopic Dirac operator spectrum entirely in terms of the 2-point function. 
The sum over topological charges of the gauge fields can be explicitly performed for these k-point 
correlation functions. A connection to an integrable KP hierarchy, for which the finite- volume 
partition function is a r-function, is pointed out. Relations between the effective partition functions 
for these theories in 3 and 4 dimensions are derived. We also compute analytically, and entirely 
from finite-volume partition functions, the microscopic spectral density of the Dirac operator in 
SU{Nc) gauge theories coupled to quenched fermions in the adjoint representation. The result 
coincides exactly with earlier results based on Random Matrix Theory. 
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1 Introduction 



While it has been known for some years that Random Matrix Theory provides universahty classes 
that describe the microscopic spectrum of the Dirac operator in theories with spontaneous breaking 
of chiral symmetry [l|, pi, [^ , it is only more recently that the precise relationship is become unraveled 
1^, |5|, P]. The essential point is that the Random Matrix Theory partition function, in a particular 
scaling limit, becomes identical to the effective field theory partition function in an analogous scaling 
limit ||l|, 0. This scaling limit is best thought of in finite-size scaling terms, as it requires sending the 
space-time volume V to infinity with the constraint that V <C 1/m^, where rrin generically indicates 
the pseudo-Goldstone masses. At the level of fermion masses rrii, one keeps the product /ij = rriiTiV 
fixed, S being the infinite-volume chiral condensate. Using the universality proof of ref. Q one 
easily establishes that the relation between the partition functions holds universally, independently of 
the chosen Random Matrix Theory potential |8|. The established identity between the two partition 
functions is not sufficient to establish that the microscopic Dirac operator spectrum can be computed 
in Random Matrix Theory, but a series of surprisingly simple relations that express all microscopic 
spectral correlators of on Random Matrix Theory in terms of the universal partition functions [Q 
indicate that one must be able to compute the microscopic Dirac operator spectrum directly from an 
effective field theory that is suitably extended by additional fermionic species. The super symmetric 
formulation of partially quenched effective lagrangians provides an analytical framework where this 
can be established |^, ^. This is an important point, because it proves, starting directly from the 
effective field theory partition function that the microscopic Dirac operator spectrum coincides with 
that obtained universally from Random Matrix Theory. 

In the derivation of microscopic spectral correlators of the Dirac operator one makes efficient use of 
the fact that one can insert factors of unity inside the field theory path integral by means of cancelling 
fermionic and bosonic degrees of freedom. Thus, for the one-point spectral function, the spectral 
density itself, it suffices to insert one such pair, and hence consider a partially quenched field theory 
partition function of the form 

Here u is the topological charge, and the additional ("quenched") fermion-boson species have masses 
m^,m^ that eventually are taken to be equal (which makes the two determinants cancel). It is also 
clear ||5| that in order to derive higher /c-point spectral correlation functions one will have to insert k 
such additional factors, which pairwise cancel in the end: 

z. = fn "./l (ti "4] ' fm. n ';!^":l! 3 ^'^^w - ^,) e--'^ . (o) 



Let us for clarity here briefiy restrict ourselves to the universality class of the chiral Unitary Ensemble 
(chUE) in Random Matrix Theory language, which is the appropriate universality class for SU(A'^e ^ 3) 
gauge groups with A'^^ fermions in the fundamental representation (as was implicitly assumed when 
we wrote down the partially quenched partition functions above). To obtain the spectral correlation 
function from the effective finite- volume partition function one needs to take a discontinuity at a cut 
in a fcth order chiral susceptibility. Using the technique of the first of ref. P] it should be possible to 
rewrite this, at the level of the effective lagrangian, in terms of a partition function extended with 2k 



additional species, all of purely imaginary masses. While this has not yet been explicitly established 
beyond the one-point function, it seems beyond any doubt that it will be possible to carry such a 
program through. The reason one can say this with such confidence is that it is known, if one accepts 
the use of Random Matrix Theory for all spectral correlators, that the A:-point spectral function in 
the chUE universality class can be be written Q 

pr (6, ....&; M) = df' n 161 n k? + ^^^ n«i - &?- — "ll"';;-''^^" . 

i = l \ /=1 / j<l Zy (l/ij) 

(1.3) 
where on the right hand side each additional mass i^j is two-fold degenerate. Indeed, the supersym- 
metric coset needed to derive the /c-point spectral function is Gl{Nj^k\k) [P, which precisely involves 
a finite-volume partition function of fc + A; additional species. 

However, it is known from Random Matrix Theory that there exists a different, and far more compact, 
expression for the spectral A;-point function: 

pJr^(6,---,a;M) = det^'^^cea,^;!/^}) , (i.4) 



a, 6 



'H/ 



where the microscopic kernel Kg {S,a,^b'i {/"«}) a-lso can be expressed in terms of finite- volume partition 
functions alone [0]: 



The formula ( |1.4| ) is just one out of three compact expressions for the fc-point spectral correlators of 
all three classical matrix ensembles, which we generically (although their history date much further 
back) shall denote Mahoux-Mehta relations Q. Taken together with the expressions (1.3) and (1.5) 



it implies a surprising identity for the partition function, which was noted in the last of reference [^ , 
and which we shall denote Consistency Condition I (to be stated in precise form in section 2 below). 

Although Consistency Condition I is valid without any doubt, it has been derived through the rather 
tortuous route of going through a Random Matrix Theory representation for the partition function. 
It is of interest to see if this identity can be proven directly from the effective partition function itself, 
without recourse to Random Matrix Theory. In this way one can logically replace the expression for 
the /c-point spectral correlation function derived through the supersymmetric technique ( |1.3| ) by the 
much more compact Master Formula (|1.5| ) and the relation (|1.4| ), without having to make use the 
Random Matrix Theory formulation. One of the purposes of this paper is to provide such a direct 
algebraic proof. 

Random Matrix Theory implies a number of other identities among the effective finite- volume partition 
functions, that cannot easily be guessed from these partition functions themselves. One of these, which 
we shall denote Consistency Condition II below, can be derived in the Random Matrix Theory language 
from the relationship between orthogonal polynomials and the kernel. We shall prove this relation, too. 
Curiously, if read conversely this one single relation allows one to derive unambiguously the effective 
partition function for any number of flavors Nf, starting with just two "boundary conditions", such 
as the effective partition function for Nf = and Nf = 1 (which are both trivial). 

Yet other partition function identities arise from the rather simple relationship between orthogonal 
polynomials in the chUE and UE Random Matrix Theories, relevant for QCD-like theories in 4 and 



3 space-time dimensions, respectively. As we shall show in this paper (see section 4), these rela- 
tions imply surprising identities among the well-known Harish-Chandra integral and the external field 
problem, both for gauge groups U{Nf). 

The existence of a long list of "miraculous" identities involving the effective partition function that 
is relevant for the microscopic Dirac operator spectrum is undoubtedly related to the fact that this 
partition function can be written as a r-function of the integrable KP hierarchy, as we shall discuss 
in section 5. The connection between the microscopic Dirac operator spectrum and this integrable 
system is certainly interesting in its own right, and may in addition be used to shed new light on the 
universal analytical expressions that have been obtained. It is also intriguing that this connection 
suggests that the effective theories relevant for describing the microscopic part of the Dirac operator 
spectrum may be "topological" in the sense of Witten ^^, i.e. having an entirely different formulation 
in terms of a BRST-exact field theory action. Intuitively this is perhaps understandable from the fact 
that the effective partition function in this regime stems from the zero modes of the pseudo-Goldstone 
bosons alone, with no kinetic energy contribution. No degrees of freedom can therefore propagate in 
this effective theory. 

Having established what we call Consistency Condition I, we get as a by-product the i^-dependent 
normalization factor in front of the spectral fe-point function. With this factor explicitly known, one 
can then perform the sum over topological charges ly of the gauge field configurations to provide a 
compact expression for this A;-point function in terms of the full effective partition functions that are 
already summed over topological charges. We do this simple derivation of the full A;-point function in 
section 6. In section 7 we turn to one of the other universality classes, labeled chSE, which is the one 
relevant for SU(A'^c ^ 2) gauge theories with Nf fermions in the adjoint representation. The analogous 
formula for the microscopic spectral density, as derived through Random Matrix Theory ^, involves 
the partition function of four additional species. We check that the resulting formula is correct at 
least in the fully quenched case by explicitly evaluating the relevant effective partition function for 
four flavors. Finally, section 8 contains our brief conclusions. Some technical details are relegated to 
the appendices. 

Before starting with the main part of the paper let us fix the notation for clarity. The finite-volume 



partition function for QCD-like theories in even space-time dimensions can be written [11| 

Zf^^d^}) = detA{M}/A{{^'}) (1.6) 

where the matrix A is defined either by 

M{f^})ij = f^i'^Iu+j-iifJ-i) , i,j = l,...,Nf , (1.7) 

or alternatively (by making use of standard Bessel function identities, and invariance properties of the 
determinant), 

Aiifi}),, = ^^l-'Ii^-'H^^^) , i,j = l,...,Nf. (1.8) 

The denominator is given by the Vandermonde determinant of the squared masses 

Nf 

Ad^n) ^ n(^'-/^|) = det [(A.2)i-il . (1.9) 

In the literature a different sign is very often chosen inside the product, which leads to an overall 
iVj-dependent sign compared with the determinant detij[{fj,jy~^] = {—'^)^^^~^'''^Tli^j{fJ-'i — ^?)- 



This just gives an overall sign in the partition function eq. (|l.6| ), which normally is irrelevant. Let 
us stress, however, that only with the definition above the partition function Zu ^ ({a*}) is a positive 
quantity. We choose the form (|1.9|) in what follows. 



2 Consistency Condition I 



First of all let us recall the consistency condition relating the determinant of partition functions with 
Nj + 2 flavors Zy ^ and a partition functions with N^ + Ik flavors as it has been stated in Q. 
Using properties of determinants eq. (11) of the third ref. Q is equivalent to 



det 

l<a,6<fc 



z'J'^iiA) 



Vi[e.-i] 



2 ^f^+''=)(M,ei,ei,...,efc,&) 



i<3 



z'J'^iiA) 



(2.1) 



The proportionality constant that remained undetermined in |Q] is therefore flxed[| as we will see below 



(fc) 



(_l)fc(^+[Af//2]), 



to be C 



In order to prove eq. ( p.l[ ) we make a more general statement which will be more easy to prove due 
to the lack of degeneracy of the fermion masses, which is present on the right hand side of eq. (p.lj): 



det 

\<a,b<h 



^^'^'^'^({/^},ea,%)' 



ri^s) 



(M) 



Vi^e. - eM - ^]\ 



^r^^^'^({M},6, 



,Ck,vi, 



,m) 



Kj 



A^f) 



(M) 



■ (2.2) 



Taking the degeneracy limit rji = ^i, i = 1, . . . , k , we recover the original claim eq. (^]^) . 



We will not be able to give a proof of the above statement in full generality. In a flrst step we 
will prove that it holds in the asymptotic regime where A*/,^j,'?/j — > cxd. In particular this flxes the 
proportionality constant in eq. (^^ which may depend on Nf and v. In a second step we can prove 
eq. (2^) for finite arguments and any /c G N in the quenched case Nf = in an arbitrary topological 
sector u. Using the established fiavor-topology duality [^ in this finite-volume scaling regime, we 
have then automatically also proven the identity for any number of massless flavors Nf in a sector of 
any topological charge u. The statement eq. (|2.2|) is this case reads 



Theorem - Consistency Condition I (massless): 

Let zi^^\{n}) = detA({//})/A({/i2}) be defined as in eqs.(|L6|)-(|r|). For Nj massless fiavors the 



following identity holds, where we have chosen the constant Zu ^ (0) to be unity: 



det 

l<a,b<fc 



zl'^^-^'k^.,,,) 



i<j 



vf)z, 



{Nf+2k) 



{Ci,---,Ck,m,---,Vk) ■ (2.3) 



Proof: (i) asymptotic region (massive, eq. ( |2.2| )): 

^ There is a misprint in eq. (11) of the third paper in ref. M], where the overall factor should read Cj '(— 1)'"(''+[^/ "!•'. 



Inserting the explicit form of the partition function into eq. ( p. 2]) and factorizing out common parts 
of the Vandermonde determinants eq. (|2.2|) can be brought into the following form: 



det 

l<a,fe<A; 



1 det A{{n},Ca,r]b) 
ill-ea detA(M) 



1 



-1-1 friS _ t2\ 
a,b=l ^ 'o- ^b) 



det A{{fi}, ^i,...,^k,Vi,---, Vk) 
det A{{fi}) 



(2.4) 



We can now apply the asymptotics of Bessel functions 



lim In{x) 



V2 



which leads to 



det^(M,U},M) 



TTX 



k (.ii+m ^s 



:(l + 0(x-i)), 



n 



n 



e^/ 



A(M,{e},M) 



(2.5) 



(2.6) 



,1 \/2vrCi??j f^-^ V'Svr/iy 
where we have taken all arguments to infinity and where A({;u}, {^}, {r/}) is now the Vandermonde 



of the unsquared sets of variables. Inserting the result eq. ( |2.6| ) into eq. (2^) we obtain for the left 
hand side (l.h.s.): 



l.h.s. 



det 

l<a,b<fc 



1 



aia+r]b 



M{fJ'},^a,Vb) 



^ (Jir+n. 



vi-eaV^^^UTb A(M) 

Nf k 



Svssnn/^'-^"'''*' -''"'./«?.» 



1 



m + ia. 



(2.7) 



and for the right hand side (r.h.s.) 



r.h.s. 



k 

n 



n 



oii+'n^ 



A(M,{e},w) 



A^l-^DtA^^f^^^i ^(M) 



^ eir+m ^f fc 



n 



n Hif^f - ^a){^^f - Va 



U.a<b(.Ca - Cb){Va - Vb) 



1 V^^^^i/JlaiV^' ""'''^■' ■'"' Ua,b=liVa + ^b 



Putting both sides together and dropping common factors we obtain 



(2.^ 



det 

l<a,6<fc 



1 



Vb+^a 



Ua<biL-Cb){r]a-7]k) 
Utb=liVa + Cb) 



(2.9) 



which is nothing else than Cauchy's Lemma. The asymptotic analysis performed so far determines 
the mass independent overall proportionality constant in eq. (^]^) to be unity. 



(ii) proof for Nf massless fermions with arbitrary u (eq. ( |2.3D ): 

Due to the fiavor-topology duality it is sufficient to prove the statement for Nf = with arbitrary 
1^ ^ 0, and then shifting u — > Nf + ly. Let us first give an outline of the proof. We will proceed with 
the simplified version eq. ( |2.4| ) of the theorem for Nf = 0, where the Vandermonde determinants 
have been already cancelled. In a first step we shall further simplify det A({,^}, {r/}) using Lemma 
1 in the Appendix W. There it is shown that detyl({^}, {r/}) is given by a determinant similar to a 
Vandermonde containing powers of ^j and 7]i as well as first derivatives with respect to these variables 



acting on a product of Bessel functions of the same index i/. In this form we can prove the theorem 
by performing a Laplace expansion of the right hand side of eq. (2.4) and using Lemma 2 for the left 
hand side. 



Starting with the right hand side the determinant can be rewritten using Lemma 1: 

k 



detA(^i,...,^fc,r/i,...,?7fc) 



fe(fe-i) 



B{0 cio 

B{v) Civ) 



n^-te)^-(^^ 



(2.10) 



i=l 



where -B(^) is the matrix of the Vandermonde determinant of squared arguments B{^)ij = {CiV ^ ^^'^ 
C(^) is the same as B with an additional ^id^^ in each row, which is defined to act only on the Bessel 
functions. If we perform a Laplace expansion with respect to the first k columns into k x k blocks we 
can use the fact that the resulting determinants det B det C can all be rewritten as the product of two 
Vandermonde determinants times k linear differential operators that can be taken out of det C. We 
obtain 



r.h.s. 



fe(fe-i) ,__. 

(-1)^ n 



1 



11 („2 _ c2\ L. 

a,6=l V'/a S6 7 o- 



^(-irA(rr^.),...,x^ 



X A(x! 



2fc 



(1)' 
2fc 



j{k)> 



cr(A;+l)' • • • ^■^cj{2k)i 



n ^-w^-.wii^' 



u\Xj 



i=k+l 



J=l 



(2.11) 
(2.12) 



where we have renamed 

Xi = Ci , Xi+k = 'ili for i = l,..,k. 

The permutations a run over all possible ( ^ ] permutations to put 2k variables into 2 sets of k 

variables with ordered indices respecting cr(l) < . . . < cr{k) and a{k + 1) < . . . < a{2k). The sign of 
the permutation is defined by (-1)'^ = (_i)i+-+'=+'^(i)+-+'^{fc) . 



The building blocks of the left hand side of eq. ( |2.4D are the 2x2 determinants 

det A{^a,rib) = [Vbdr,, - iad^.)h{ia)Iu{rib) (2.13) 

where we have used again Lemma 1 eq. ( [A.2| ). Therefore we can rewrite the left hand side as 



l.h.s. 



det 

l<a,b<k 



rit - a 



det A{^a,Vb) 



det 

l<a,6<fc 



— ^ {Vbdrj, - Cad^ 



vt - a 



n^-teK-(^0, (2.14) 



J i=l 



with the derivatives only acting on the Bessel functions. 



Using again the fact that determinants differing only by a single column can be added, eq. ( |2.14 ) can 
be rewritten as a sum of determinants containing only one derivative %Q^j, or ^ad^a^ ^ ~ ^, ■ ■ ■ ,k. Due 
to the structure of the determinant in the columns with ijbdr^^ the derivative can be taken out as a 
common factor. After reordering columns we end up with 



det 

l<a,6<fc 



Ik k 

^(%4-^a4) = EE(-ir'^^detG(.;7.) U ^7.a)^..,u)' (2-15) 



vt - ea 



i=0 1i 



j=i+l 



where the matrices G{i;'^i) are defined in Lemma 2 eq. (A. 11) to contain derivatives £,ad^^ in the first 
i columns. The 7^ are all ( ^ 1 permutations for different G{i;ji) with i fixed, where the ordering is 



such that 7i(l) < . . . < 7j(i) and 7j(i + !)<...< 7i(A;). The factor (—1)' stems from taking out the 
common factor minus one of all columns containing derivatives with respect to the ^'s. 

We can now apply Lemma 2 where det G{i; 7i) is evaluated by Laplace expansion and Cauchy's Lemma, 
to obtain 



det 

l<a,6<fc 



vt - e 

k 
fe(fc-i) .__. 

(-1)^ n 



3T2 (%^% - ^^^ia 



11 („^~W) 5m(~i)'^'^'^'^*^(^?»(i)'---'^?,(fc-i)'^7«(i)'---'^iw) ^ 

a,b=l ^ 'o- ^b) j=0 7i.7i 

k k 

X M^i(k-i+l)^- ■ ■ > 4(fc)' ^7.{i+l)' • • • ' ^7»(fc)) n ^7.(J)%,0) n '?7.(0^^7.(0 (2.16) 



This equation is nothing else than the operator in eq. ( 2.11| ) acting on nj=i Iu{^i)Iu{Vi)> where here 
in the first sum the number of ^'s in the first Vandermonde is made explicit to be i. 

To see this we observe that in eq. ( p. 161 ) we have X]j=o { i ) ~ { k ) different terms, which matches 
to the number of permutations in eq. ( |2.11 ). In order to map individual permutations including signs 
we use again the notation eq. ( 2.12[) . The index 7(j) of the r/'s thus changes according to 



7(j) -^ k + jij) = a{j + k) , 
which implies(-l)T(i)+-+'^W = (_i)<T(i+A:)+...+<7{i+fc)-ifc^ Consequently we obtain 

(_l)7i + li + i — (_i 



(2.17) 



^^^^ + Et">«(fc-o 'ik+ "'-"<,"-'+^> +Et";7.(0 + * 



fc(fc+l) , V^fe , .N 



(2.18) 



where cjj is a permutation with i ^'s with ordered indices {7i(j) = cri{j); j = 1, . . . , k — i} £ {1, . . . , k} 
and (k — i) ry's with ordered indices {(Ji{k+j); j = l, . . . ,i} G {k+1, . . . , 2k} as in the first Vandermonde 
in eq. ( |2.1l| ) . We have thus completed the matching of both sides of theorem eq. (^^) . 



3 Consistency Condition II 

In this section we will prove the following theorem relating partition functions with Nf, Nj + 1 and 
Nf + 2 massive flavors, as it has been stated in the last of reference Q. 

Theorem - Consistency Condition II: 

For zi^^'\{fx}) = det A({/i})/A({^2|) ^g defined in eqs.(|]|)-(|rD it holds 

^r^'^^^({M},e,.) 



(e2-r?2)^f/)({^}) 

Nf 






J2^^^^,^+C^^)zl'''^'\{^,},0\zi'''^'\{^^},rl) - (e-^) 



(3.1) 



A few remarks can be made here. Taking the inverse statement, the differential equation eq. (3.1) 
together with the boundary conditions Zu = 1 and Z{, (/i) = /i/(/i) can be seen as a generating 



equation for all Zy ^ ({/"})• In this way we can actually derive the precise form of the partition 



function eqs. ( |1.6| )-( [L9D instead of taking it as a starting point. In ref. |13|, an even more compact 
recursive relation for the partition functions was derived in the Random Matrix Theory formulation, 
using the supersymmetric formalism. 



The consistency condition eq. ( |3.1D had been conjectured in the third reference of [Q] (eq. (16)) with 
an arbitrary constant, which has been determined here to be C = 1. This result follows independently 
from inserting the asymptotics of the Bessel functions. 



Proof: Before going into the details let us give the general outline of the proof. We will first investigate 
the action of the power-counting operator Yl,il\ P^i^^in +Cf^C on the partition function Zy ^ ({/u},^). 
Most of the outcome will be again proportional to the same partition function which then drops out 
due to the antisymmetry with respect to ^ and r/ in the bracket in eq. ( |3.1D . The remainder times the 

partition functions Zv ^ ({A'}) '^) expanded once will then precisely arrange to the Laplace expansion 
of the left hand side with respect to the last two columns. Throughout the proof we will make use of 
3 different Lemmas collected in Appendix P. 



Let us start with the action of the power-counting differential operator on the Vandermonde in the 
denominator. It is easy to show that 



Nf 



(E^^^M.+e95)A({^2|^^2) ^ (Ar^ + i)Ar^A({^2|^^2) 



(3.2) 



i=\ 



since the operator counts the sum of all powers. From the product rule every factor in A gets 
differentiated twice and is thus reproduced with a factor of 2 in front, furthermore there are [N j -|- 
l)Nf/2 such factors. Hence the differentiation of the A's drops out of eq. ( |3.1D due to the (^ <-^ rfj 
antisymmetry. Inserting the explicit form of the partition function we obtain 



det^(^,r/) = -(detA)-^ 



i=l 



(E /^i^M. + idi:) det A{i) det ^(r/) - {^ ^ r,) 



(3.3) 



after cancelling the Vandermonde determinants. Here and in the rest of this section we have omitted 
the dependence of the matrix A on the set of variables /Xj. Next we apply the differential operator to 
the determinant of A{^). Since it is linear and the matrix A depends on each variable only in one row 
we can differentiate row-wise, using 



^JL^^[^J^Iy+n{^^)) = (2n + i/)/x"I^+„(;u) + /i"+^I^+„+i(;u). 



(3.4) 



In fact this is the only place where the properties of Bessel functions enter, the rest of the argument 
being valid for general matrices A. Determinants that only differ by one row can be added and we 
thus have 



E 



IXid^- det A 



Nf 



k=l 



Nf 



2 ^ det Afc + 1/ det A + det Ak 

L=2 



(3.5) 



Here in the last term Ak is the matrix A with the k-th row shifted to the left by one unit, {Ak)ki = 
fi^Iy+iink) and Ak = A else, which results from the last term in eq. (^). This last term which can 
be further simplified using Lemma 3 in Appendix ^ will be the only term that survives in eq. (| 
The second term u det A vanishes immediately due to the antisymmetry in eq. ( |3.3| ) . 



We will now show that also the first term in eq. ( |3.5D is proportional to det A and thus drops out. The 
matrix Aj^ is defined to be the matrix A with the first L entries in the fc-th row vanishing, {Aj^)kin = 
for m = 1, . . . ,L and Aj^ = A else. The sum over L in eq. (|3.5D thus reproduces the 2n//"/jy+„ (^) 
from eq. ( |3.4[ ). Applying Lemma 4 from Appendix^ we obtain 



Nf Nf 



N 



f 



2^^det^^ = 2^(iV/-L)det^ = {Nf - l){Nf - 2) det A . 

L=2k=l L=2 

Together with Lemma 3 Appendix ^ we finally obtain 



(3.6) 



4 = 1 



Nf Nf 

J2f^idf.,detA = ((iV/-l)(iV/-2) + i.Ar^)detA + X](-l)'^^+Vf'^.+iV,(Mj)detA*^^. (3.7) 



i=i 



where A*j^ is the matrix A with row j and column Nf missing (algebraic complement of matrix 
element Ajn ). Inserting eq. ( p.7| ) for the derivative of the (A'^^ + 1) x {Nf + 1) determinant of A{^) 
in eq. ( |3.3D we obtain 



Nf 



det^(C,ry)detA = (^ ^(-l)^/+i+Vf''^'/.+JV,+i(/^i) det A*^^.+i(0 



- e'^/+i/,+iV,+i(6det^jdetA(r?) - {^ ^ t]) . (3.8) 

Expanding det A{r]) and det A{S,) on the right hand side with respect to the last column and multiplying 



out we obtain for the right hand side of eq. (3.8) 



Nf 



Nf 



r.h.s. = J2 (-l)*^Vj ^ Iy+Nf+i{lJ'j)lJ'i ^ly+Nf (Mi) [det ^*jv^+i(0 det A*^^+i(7y) - (^ ^ ry) 



Nf 



+ [ Y.{-lff-'^^^li'^'^'h+Nf+i{H) (^^^■4+7V,(^)detA*^^+i(0 -e'^^/.+iV,(e)detA*^^.+i(r? 

Nf 

^(_l)iV/+i+i^^^//^^^^(^^.)(^A^/+ij^^^^^^(^)det^*^/r?)-ry^/+i/,+;v,+i(^)detA*^/0 



i=i 



+ 



[e^^^Iu+Nf+mv'''Iu+Nf{v) - r?^/+i/,+^,+i(r?)C^//,+^^.(0)detA 



det A 



(3.9) 



We are now ready to apply Lemma 5 from Appendix ^ which relates products of determinants of 
matrices which just differ by the last two rows. The bracket in the first line then reads 



det^*^^,+i^,.^^.+2(e,r/)det^| ;j+|| ^ < ^^ 



(3.10) 



det ^*;v,+i(6 det ^*;v,+i(^) - {i ^ v) 

where we have introduced the matrix A*jq _|_i.,-^ _|_2 where i-ih. and j'-th row and the last two columns 
are missing. The obvious symmetry A^jq _|_i.,-7Vf+2 — -^*jN ■+i-iNf+2 leads to a sum over i < j only in 
the first line of eq. (|3.9D . Since we can trivially rewrite A*j^j^-^{£) = ^*7v +i-JV +2Af +2(?'^) °^e '^^'^ 
immediately see that eq. ( ^.g| ) is nothing else than the Laplace expansion of det ^(,^, r/) in 2 x 2 times 
Nf X Nf blocks choosing the last two columns for the expansion. Thus we have obtained the right 
hand side of eq. (^^) and completed the proof. 
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4 Relations between QCD3 and QCD4 partition functions 

Let us first recall how the Consistency Condition I and II of the the previous sections have been 
originally derived in Q. The important point is that every finite volume partition function equals a 
partition function of a massive chiral random matrix model in the microscopic large- A^ scaling limit. 
The fact that in matrix models the orthogonal polynomials, the associated kernel, as well as all cor- 
relation functions can be expressed in terms of matrix model partition functions translates the well 
known relations between these objects to relations among finite volume partition functions. Consis- 
tency Condition II reflects the Christoffel-Darboux identity between the kernel and the orthogonal 
polynomials, and Consistency Condition I is the Mahoux-Mehta relation among fc-point correlation 
functions and the determinant of the kernel. 

In this section we will follow the same reasoning to translate relations between matrix model quantities 
for QCD3 and QCD4 to the corresponding finite volume partition functions. In contrast to the previous 
section we will not provide a proof solely based on finite volume partition functions. While the 
identities are very easily derived in terms of the Random matrix Theory formulation, they translate, 
in the microscopic large-A^ limit into highly non-trivial relations between group theory integrals of 
Harish-Chandra type (for unitary groups) and what can be called the external field problem for unitary 
groups. These surprising identities deserve to be understood in their own right. 

Let us start with a relation between the orthogonal polynomials the chUE (QCD4) and those of the 
UE (QCD3): 

P\.%u7""\-'-M) = P^nAe M^}) • (4.1) 

This simple identity is very easily derived from the defintion of the two Random Matrix Theories 



as given in [|T6|. Since in ref. [l^ no explicit use was made of the measure, eq. (4.1) also holds in 
the massive caseg. On the right hand side the same Nf masses appear in pairs with opposite sign. 
Expressing the orthogonal polynomials by partition functions as given in the second of ref. [Q] we 
obtain 

''"^^ e'i/.(M) 4S({M)) ■ *■* 

Here the unknown proportionality constant only reflects the choice of normalization for the polynomi- 
als, and it can easily be fixed. Since we have started with an even polynomial the odd-flavor partition 
function on the right hand side is given by Q 

4S'^ = j dU co^h[TT{MUTU^)] , (4.3) 

where 7W=diag (/xi, . . . ,//Ar^, — /^i, • • • , —fiNf,iO ^^^ r=diag(lAr^. ,-lAr^,_l_i). The denominator is given 
by the 3-dimensional even-flavor partition function [^, 

^qSdI = JdUexp^IiiMUrU^)] , (4.4) 

with A^=diag (/ii, . . . ,fJ-Nf, — /^ij ■ ■ ■ , —fJ-Nf) and r=diag(l jy^ r'^Nf ) ■ These group integrals are both 
of the Harish-Chandra type, while the left hand side of eq. (|4.2| ) is given by an entirely different ratio 
of unitary group theory integrals of the external-field kind. This is the first of such relations. 

■^In contrast to llq] we have shifted the non-integer part to the left, u = —1/2, in order to deal with a physical (even) 
number of massive flavors for QCD3 on the right hand side. 
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The relation between the orthogonal polynomials eq. (^]^) can be exploited furthermore in order to 
relate also the corresponding kernels of the chUE and the UE. For convenience, let us introduce the 
"wave functions" 

^„(A) = ,/Z{X)Pn{X) , (4.5) 

where uj{X) is the measure function (so that the wave functions ^n(A) are orthogonal with respect to 
a weight of unity). Now, from eq. (3.11) of ref. |]l6[ we can use the wave functions of the UE inside 
the Christoffel-Darboux identity for the chUE kernel^. Doing this, we readily derive the following 
identities: 

^N, chUE [^ ^W ) — 2 _ y,2 \^2N{^) ^ ^2N-l{^) ~ ^ ^ 2N-l{^)^2N {^) 

1 C2Ar (ir,UE,Mv,UE /n ,j,UE i^mt.UE 



^ ^2JV (,^,UE(^iv,UE / N ^UE (^\^UE 



[W, 



5 (tI^ (*"" <"'*^"-'<"'> + *'^«-i(^)*'^«(» 



\{- 



(2^/) (^ .„^ , I<'(27V;) 



(z,w) + Kf^fljEi-^^^)) , (4.6) 



2N, UE 

where for clarity we have not explcitly indicated the mass dependence of the wave functions. The coef- 
ficients C2N are defined in ref. ||l^ . We have here made use of the fact that in the UE the polynomials 
P„(A) are of parity (— 1)*^. As one can easily convince oneself the last equation is, despite its appear- 
ance, symmetric in the arguments z and w. Expressing the kernels in terms of the corresponding finite 
volume partition functions 0| eq. (|4.6D leads to the following relation: 



(—1)' 2 1 2 y'^LU _ I -^ -^ 



(4.7) 
The proportionality constant is known in this case as it can be determined for both kernels from the 
matching condition with the microscopic spectral density. 

Let us finally mention that also in the case of QCDs-like theories with an odd number of massive fiavors 
a relation similar to eq. ( [4.6| ) exists. In ref.|14| the random matrix model kernel for odd-flavored QDC3 



has been derived from a chUE, which in the microscopic scaling limit reads: 



K'siHUE'i^^^-A^}) = K'sJEit^-Af^},0) . (4.8) 



Inserting again the representation in terms of partition functions |Q], we obtain the following relation: 

^ ' ^Zijsp)) " ■"- «r'(w.o) ■ 

For more details on the odd-flavor partition function we refer to ref. p!^ . 

We have explicitly checked all of the above relations for a few (small) number of flavors, starting 
from the finite-volume partition functions alone. Apart from those given above, we have empirically 
found other non-linear relations between these Harish-Chandra type integrals and those of the unitary 
external field problem. These relations, however, do not seem to follow easily from the Random Matrix 
Theory formulation. 

■^Eq. ( |4.l| ) trivially holds also for the wave functions instead of the polynomials. 
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5 The effective QCD4 partition function as a r-function 



The identities derived above may have their origin in a surprising relation to integrable systems, which 
we shall now briefly discuss. The starting point is the following observation. Suppose we deflne an 
Nf X Nf hermitian matrix integral by 



t{X) = f dYexp[TiiXY + V{Y)]] . 



(5.1) 



With X itself being an A'^^ x Nf hermitian matrix, and the potential V{Y) as yet unspecified, this is 
a generalized "external field problem" of Random Matrix Theory. The universality of its correlation 
functions have been proved in the microscopic large- A^j limit in [O] for polynomial potentials V(Y). 
Surprisingly, a closed solution of the integral eq. ( |5.1[ ) can be written down for any value of Nj, and 



for any potential V{Y) that satisfies suitable convergence criteria |17, 18|. After diagonalizing the 
y-matrix, Y = [/diag(yi, . . . ,yNf)U'^ , one obtains the standard Jacobian of A(y)^, and one can then 
make use of the Harish-Chandra integral to obtain 



r{X) 






J2{x,y, + Viyj)) 



(5.2) 



where the Xi are the Nf eigenvalues of X. If one now introduces the function 



as well as the derivatives 

one sees that the integral is simply 



4>{x) 

Qk 



dx'^ 



t{X) 



_ de,i[4)j^i{xi)] 
A(x) 



(5.3) 
(5.4) 

(5.5) 



This shows that t{X) is a r-function of the integrable KP hierarchy. 

The expression ( |5.5| ) has an uncanny resemblance to the finite- volume partition function if one identifies 



Xi = fif. This is particularly clear if one considers the form (1/7), and starts with the case u = 0. 
Using the Bessel function identity 



d" 



(x^//^,(x)) = -L^^/-fci^^_,(2;) 



we see that the partition function (|1.6[) can be written in the form (|5.5|) if we identify 



Hf^') 



Nf 



^^ ^NfiytJ''^) , 



(5.6) 



(5.7) 



and ignore irrelevant overall factors. The case of non-zero ly can then be treated by using again the 
flavor-topology duality |I^, thus obtaining the v ^ case by simply setting v of the A'^^ masses to 



zero. In fact, eq. (5.7) alone suffices to prove that the effective partition function is a r-function of 



the integrable KP hierarchy |17]. It is nevertheless interesting to note that in addition an integral 
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representation actually exists such that the partition function explicitly can be written in the form 
( ^.1| ) [|l9|. This turns out to correspond to a potential V{Y) = 1/Y — A'^^ln(y), and an integration 
contour for the eigenvalues encircling the origin. (Strictly speaking this is outside the scope of the 
herniitian matrix formulation (5.1), so the notion of "hermiticity" is here simply taken to mean "of 



flat measure" - see ref. [19| for a discussion of this point). 



Knowing that the partition function is a r-function immediately implies a number of identities. The 



most general of these is the following set of Hirota equations, which read |17| 



= {Xa-Xb)T{X;Pa,Pb,Pc + 'i-)T{X;Pa + l,Pb + 'i-,Pc) 
+ {Xb - Xc)T{X;Pa + l,Pb,Pc)T{X;Pa,Pb + l,Pc + 1) 
+ (Xc - Xa)T{X;Pa,Pb + l,Pc)T{X;Pa + l,Pb,Pc + 1) 



(5.8) 



Here the pi denote the multiplicities of the parameters Xi for i = a,b, c, where in our notation Xi = fJ-f- 
Due to the Jacobi identity for determinants yet another identity holds for r-functions as given in e.g. 



eq. (2.43) ofref.[ll7[: 

T-(^/+2) (-{3.}^ XNf+l,XNf 



{xNf+1 - XNf+2)T^^f\{x}) 

r(^/+i)({x},x^,+i)f(^/+i)({x},x^,+2) - (X^,+1 



XNf+2) 



(5.9) 



The upper index indicates the number of parameters of the corresponding r-function and f means that 
in the last row the index of the functions (pkix) in eq. (5^) has been shifted by +1. This relation looks 



remarkably similar to our Consistency Condition II eq. (3.1). However, in eq. ( |3.lD the derivatives, 
which shift the indices of the (pk are taken with respect to all variables. In the derivation of the 
Consistency Condition II we have used properties of the Bessel-functions in only one step, namely in 



eq. (3.4). The fact that all r-functions eq. (5.5) obey the property 

2/i5^0fc(a;) = 5x(/>fc(2;)U=^2 = (pk+iix) 
can probably be used to show that all r-functions obey our Consistency Condition II. 



(5.10) 



There exists another set of relations for r-functions of the KP hierarchy which apparently can be 
related to our finite- volume partition functions. According to ref. pO] these relations read 



det 

l<a,6<fc 



r{t + ICa'] - IVb']) 
{^a - %)r(t) 



k 
a<b 



^b)iVb-r]a) 



r(t + E'a=li[Ca']-[ria']) 



UUCa-VbMt) 



(5.11) 



At first sight they look remarkably similar to our Consistency Condition I in the form of eq. (|2.2| ) . In 
order to to explain the differences let us give the notation of eq. ( |5.11| ) from ref. [20|. The argument of 
the r-function t stands for all the coupling constants or times in a matrix potential V{X) = J^'k^i ^feA . 
The bracket [ 1 then is a shorthand notation for 



r(t ± [C] 



T{h±C,t2±le,t3±k^...) ■ 



1 

3^ 



(5.12) 



Performing the sum over the additional parameter [^ ^] in the potential leads to an additinal loga- 
rithmic term V{\) -^ ^(A) — ln(l — X/S,), which resembles an extra "mass term" if we could consider 
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this as an ordinary hermitian Random Matrix Theory in which one has taken the microscopic hmit. 
Such an identification is, however, far from obvious. Moreover, in eq. ( |5.11| ) these extra terms occur 
in pairs ^ and r] with opposite signs and thus one of them appears as a bosonic "mass term". Due 
to this difference, apart from other additional factors, we have not been able to explicitly match our 
Consistency Condition I eq. (|2.2D with eq. ( 5.11| ). Moreover, no direct proof has been given in ref. 



|20[ | for the relation eq. ( 5.11| ), and we have not been able to find it elsewhere. 



6 Summing over topological charges 

So far our discussion has been restricted to finite-volume partition functions Zi, ^ ({a*}) in sectors of 
fixed topological index v. These partition functions can be thought of as Fourier coefficients of the 
full partition function, which for given vacuum angle 9, is given by 

oo 

Z(^f\e,{f^}) = E e-^zf ^) ({/.}) , (6.1) 



iy=—oo 



or, in terms of the effective partition function on the coset of chiral symmetry breaking in this case, 

Z^^f\e,{fi}) = f dUexp\v^Re[e'^/^fTrMU^]\ . (6.2) 

JsU(Nf) L J 

Contrary to the effective partition functions Z^ ^ ({^}) in sectors of fixed gauge field topology, the 
group integral of eq. (^]^) is not known in closed form for Nf > 3. We are not aware of any analogue 
of the theorems discussed above for the full partition functions, and in view of the non-linearity of the 
relations it seems unlikely that they could be established. Nevertheless, as an interesting by-product 
of the above analysis we are now able to provide a simple compact formula for any A;-point spectral 
correlation function after having summed over all topological charges. 

The first observation is that for any observable (O)^ in the fixed-z^ theory one finds the same observable 
in the full theory by summing over v with weight factor e Zp ^ ({/^}): 

oo 

{{O)) = E e-^zf ^)({M) (O). . (6.3) 

i'=—oo 

Next, noting that the /c-point spectral correlation function by itself is just an expectation value, also 
this function can be summed over topological charges: 

oo 

Psi^i,...,^k;0,M) = 2^'''\O;{f,})-' E e-^zf^\M)/>lr^(ei,...,a.;M), (6.4) 

^■= — 00 

where we have not included the zero-mode contributions in the spectral sums, and where we have 
already taken the microscopic limit. 

We now make us of the fact that the A;-point function in a sector of fixed topological charge v can 
be expressed in terms of a partition functions with 2k additional species as in eq. ([L.3|). As a 
side result of proving the Consistency Condition I in section 2 we have already fixed the constant 
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6*2 = (— 1)'^('^+[^//2]). Only the i^-dependence is important for the summation over topological 
charges, where it leads to a shift in the 0-angle: 



Nf 



P5(ei,...,e.;e,M) = z(^f\e;M)-' Y. e-^-i)'^('^+[^//2i)n led n(e.' + /^/)i X 



/=1 



3<l 

k ( Nf \ k 

(-i)'=[^^/^in le^iKe^'+M/) iK^i-ezY 

i=i V /=i / j<i 

Z(^/)(e;M) 



X 



(6.5) 



Due to the periodicity of the angle we need to know the partition function at either a shifted or 
unshifted vacuum angle 6 for the /c-point correlation function with k either even or odd Q 

7 Fermions in the adjoint representation 

With Nf fermions taken in the adjoint representation of the gauge group SU(A''c > 2) the pattern 
of spontaneous chiral symmetry breaking, if it occurs at all, is believed to proceed according to 
S\]{Nf) — > SO(A^/). In the Random Matrix Theory classification this corresponds to the chiral 
Symplectic Ensemble chSE. Fermions in the adjoint representation of the gauge group occur in e.g. 
supersymmetric gauge theories even without matter fields, but its interest in the present context stems 
more from the fact that staggered fermions in the fundamental representation and gauge group SU(2) 
actually also fall into this universality class away from the continuum limit. In ref. [22| the effective 
partition function in the same finite- volume limit as above was written, in a sector of fixed topological 
charge v = NcU, in terms of a relatively simple- looking group integral over the unitary group \]{Nf): 

zf^\M) = /"(iC/(det[/)-2^exp[sFReTr7W[/[/^] , (7.1) 

where as before M denotes the mass matrix. This group integral is surprisingly difficult to perform 
explicitly, and it has in fact until now only been evaluated for Nf = 2 fermions of equal mass, and any 
topological charge i^ [23|. In general, the equal- mass partition functions are much easier to evaluate 



due to an interesting rewriting of the group integral ( |7.1| ), which is valid for those cases |22|. We shall 
here use that form to explicitly evaluate the group integral for Nf = 4 equal-mass fermions in a sector 
of zero topological charge (the derivation extends straightforwardly to any topological charge, but we 
have not considered that extension in detail). 



The partition function ( 7J ) of Nf equal masses has been be rewritten by Smilga and Verbaarschot as 
(for Nf even) |2|: 

zff\M) = Pf(A) , (7.2) 

where, in our normalization, the Nf x Nf matrix A has elements 

A„ = -'- r ^ r ^e(0-(/<)e^(P*+5^)e'^^°^*+^™^^+*^('^+^) , (7.3) 

2 J_7r 27r J_7r 27r 



In pi| a factor of ( — 1) ^ is missing in the formula for the density (corresponding here to fc = 1). 
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and the indices p and q run from — 5^ + 2 ^° ~2' ~ 2' ^S^i^ 1-'^ — mTiV is the (common) rescaled 
mass. A for our purposes convenient infinite-sum representation was also given in ref. ||2^, based on 



a Fourier series expansion of the sign function in eq. (7.3): 



00 1 

^P9 = E ^^j^:p[h+p+k^ip)h+q-k-\{^J) ■ (7-4) 

k=—co 

It is with this form of the matrix A that we have managed to evaluate the partition function (^]^) for 
A'^^ = 4 equal masses and P = (this last restriction can readily be lifted). After some tedious algebra 
we find (technical details can be found in Appendix C): 

4'^/^) = ^ \h{2lif - h{2^^f\ + :T^/o(2/i) /'''(it/o(2/x) . (7.5) 

The last integral is explicitly known in terms of a combination of Struve and Bessel functions, but we 
leave the result like this in order to facilitate a comparison to be discussed below. 



The reason for our interest in the partition function (7.1) is a general relation derived in the third of 
ref. 0], which expresses the microscopic spectral density of the Dirac operator for this case in terms of 
the partition function itself and the partition function with 4 additional fermion species of imaginary 
(degenerate) masses: 



S^)(c i.w - n c3 (c2 , ,.2^4 z j^^^^ ({MJ^I^ g}) 



pVii; {^}) = c, e ie + ^^r '^ ^^n,; ■ (^-^^ 



The normalization coefficient C4 can be fixed as soon as one settles on the normalization of the partition 
functions. The general formula ( [7. 61) as derived in ref. Q in the Random Matrix Theory formulation 
(and the partition functions involved were therefore those of Random Matrix Theory, too). But in the 
microscopic limit these coincide, modulo uninteresting mass-independent normalization factors, with 
the field theory partition functions ( |7. 1| ) , thus giving explicitly the microscopic spectral density in 
terms of the field theory partition functions, as indicated. Although this general formula may provide 
a simple way of deriving the massive double-microscopic spectral density for this universality class, it 
has not yet been tested due to the lack of a simple analytical expression for the partition functions 
([7. 1|) . Now, with the analytical result (7.5) we can for the first time check the formula, since for 



the quenched case (formally defined by taking Nf to zero) the partition function itself becomes an 
uninteresting constant (which we take to be unity), while the partition function in the numerator of 
eq. ( |7.6D is that of just four fermions (of imaginary and degenerate masses). 

Using well-known relations between Bessel functions and modified Bessel functions, we thus derive the 
quenched microscopic spectral density of the Dirac operator for this case: 



pPio = c,ezi,'\{ic}) 



= cAaM2Cf + Jii^Cf] - ^Jom I dtJo{t)\ . (7.7) 

If we next impose the matching condition /9s(C — >■ 00) = l/vr, then the overall constant in front is fixed 
to C4 = 1. This finally gives 

pPio = ^Jomf + Jimf] - ^Jom rdtut) , (7.8) 

^ JO 
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which agrees exactly with the result obtained directly from Random Matrix Theory |24]. 



In a derivation of the formula (W^ directly from a partially quenched chiral Lagrangian, the four- 
fold mass-degenerate additional fermion species should come out from coset of the supergroup chiral 
symmetry breaking of that case. The fact that one is led to four additional quarks is a somewhat 
surprising feature, as in a very naive counting one could have expected two: one from the additional 
quenched quark, and one from its supersymmetric partner. It is therefore very comforting to see 
that the recent explicit computation of the partially quenched effective lagrangian by Toublan and 
Verbaarschot [^ in this case precisely leads to four additional species in total. By compactifying 
variables after taking the discontinuity that gives the spectral density it should therefore now be 



possible to derive the formula (7^) directly from the effective Lagrangian, following the steps of the 
first of ref. §. 

8 Conclusions 



Our main purpose here has been to show that the surprising relations among the effective partition 
functions relevant for describing the microscopic Dirac operator spectrum can be derived directly, 
without recourse to the Random Matrix Theory formulation. We have noted that these identities 
most likely have as their origin the fact that the effective partition functions of the chUE universality 
class are r-functions of an integrable KP hierarchy. We believe that one of these identities, here called 
Consistency Condition II, holds in general for all these r-functions. 

As a by-product of our analysis, we have computed the z^-dependent normalization factor of the k- 
point spectral correlation function of the same universality class. This has allowed us to perform the 
sum over topological charges ly explicitly, and express the /c-point function of the full theory entirely 
in terms of the full effective partition functions, without the restriction to fixed topological charge. 

We have noted a series of relations between the effective finite-volume partition functions for QCD3- 
like theories and QCD4-like theories. These relations translate into surprising relations between the 
external U{Nf) field problem and the Harish-Chandra integral for unitary groups. 

Finally, we have considered an analogous formula for the microscopic spectral density of the chSE 
universality class, which expresses this spectral density in terms of the effective field theory partition 
function with four additional (imaginary-mass) fermion species. We have explicitly shown that this 
formula yields the same analytical result as the Random Matrix Theory approach in the quenched 
case of Nf = 0. Considering the analytical difficulties in extending the corresponding Random Matrix 
Theory calculation to the case of massive fermions, this may provide the most economical way of 
deriving all the microscopic spectral correlators of that universality class. 
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A Some lemmas for Consistency Condition I 

Lemma 1: Let A ({^}, {r/}) be the following 2k x 2k matrix 

(^({e},w))„- = <! 'ur^'-'H . t:v'\, vj. (a.i) 



erVj_ite) i = i,...,k 
"' ^ vr'i.+j-iivi) i = k+i,...,2k 



Then the following statement holds: 

del Am, {V}) = (-l)^('^-i'/2 



B{0 cio 

B{v) Civ) 



n^-te)^-(^0 (A-2) 



i=l 



where the two k x k matrices B and C are given by 

B{O^J = iily-' , C{i)^, = {e^y~H^^., i,j = l,...,k . (A.3) 

In other words B is the matrix inside the Vandermonde with squared arguments whereas C contains 
in addition the power counting operator in each row, which acts on Bessel functions outside the 
determinant only. 

Before proving the lemma let us add a remark. The above statement is not specific for a matrix with 
variables split into 2 equal groups, as we will see in the proof below. Furthermore eq. ( [A.2| ) trivially 
extends to matrices A of odd size, where one additional column has to be added in the matrix B. We 



have just presented the statement in the precise form as we need it in the proof of theorem 2.3. 
Proof: Let us first for simplicity set 

Xi = Ci , Xi+k = m for i = l,...,k . (A.4) 

We will in a first step use the property of the Bessel functions 

Xln+l{x) + 2nln{x) = xIn-1 (A.5) 

to reduce all indices of the Bessel functions in detA({j;}) down to v oi u + 1. Starting with the last 
column we can add 2{v + 2k — 2) times the last but one column to it to reduce the index of the Bessel 
function by 2: 

xf-^I,+2k-i{xi) + 2{u + 2k-2)xf-^h+2k-2{xi) = xf-^I,+2k~z{xi) . (A.6) 

We proceed similarly with the last but one column and go down till the third column. All indices of 
Bessel functions have been reduced by 2 except in the first 2 columns. We then start again with the 
last column reducing the index by 2 (eq. ([A. 61) for z/ ^ z^ — 2) and continue down to the fifth column. 
So after going down through all columns k — 1 times always starting from the right we have achieved 



det^({x}) = ^ det x| h ,^^^^{xi)] . (A.7) 

Li a second step we write all functions Iyj^i{x) which appear in every second column as a derivative 
dxlv{x) so that we can factor out all Bessel functions to the right of the determinant. Using eq. (3.4) 
for n = we have 

x2'+il,+i(x) = x^^+^dj,{x) - ux^^Iuix) . (A.8) 
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Inserting this expression into every second column we can again compensate the last term in eq. (AJ 
by adding v times the column to the left. We thus have 



detA({x}) = deti({x})]J/^(xi) , 



where 



(^({^}))^, 



i=l 



T-C ■^T^'f'V'^l, ^ Vi = i,...,2fe. 



A. d^,. 



j = 2,4,...,2/c 



(A.9) 



(A.IO) 



Here we have pulled out all common factors luixi) out of each row i and again the derivatives only 
act on these. After reordering those columns with even powers to the left and those with odd powers 
times a derivative to the right, which results into a factor (— l)^'('^^i)/^^ we obtain eq. ([A.2|) when 



renaming back the variables from eq. ( A.4 ). 
Lemma 2: Let 



{G{^ ■■ l)h 



^e«4 



n^ — ^f^ 



l = i + l, 



.k 



V j = 1, . . . , /c 



(A.U) 



where 7 is a permutation of the indices 1, . . . ,k with the ordering 7(1) < . . . < 7(1) and j{i + 1) < 
. . . < ^{k) for 1 < i < k fixed. It then holds 



fc(fc-i) 
detG(i;7) = (-1)^ 



1 



11 ('„2_t2 
a,b=l ^ '<^ ^b 



2_^(-l)'''A(^^(^),... ,^^(fc_j),T/^(i),... ,r/^(j)) X 



X A(^^ffc_j^i^,...,^^ffc^,r/^fj^i^,... ,?7^(^)) Yl ^7(j)% 



57(fe-i+l) ) ■ ■ ■ 5 ?7(fc) ' 'l-y(i+l) ' 

where 7 is a permutation as 7 of the indices of the ^'s. 



-lU) 



(A.12) 



j=k—i+l 



Proof: We will prove Lemma 2 by doing a Laplace expansion with respect to the first i columns into 
i times {k — i) dimensional determinants. Each of the subdeterminants can then be evaluated by using 
Cauchy's Lemma eq. (|2.9|) for squared variables, which reads 



det 

l<a,fe<A: 



1 



We thus obtain for the first upper block 



(-1) 



I{i,=Ma-il) 



det 



1 



Hi) ^i 



2^19^1 



(-1) 



'('-!) Ila<biCa ~ ^b)\V^(a) ~ ''l-y(b)) 



Ila,b=li'n^(a) ^b 



7(a) '^(^)^n^.% 



(A.13) 



(A.U) 



j=i 



and similarly for the lower blocks, where we can directly use eq. ( |A.13D . Taking all the permutations 
7 according to the Laplace expansion we arrive after a few steps at 



detG(z;7) 



(_l)%^(_l)i(i-fc) Y[ 



ah=l(^a 4) 



2_^(-l)^A(^^(^), . . . , ^^(j) , T/^(j+i) , . . . ,r?: 



-y{k)j 






(A.15) 



20 



where the sign of the permutation is defined by 

(_1)7 = (_i)E;=i(i+7(i)) ^ (_i)E-=.+i(j'+^(j')) . (A.16) 

To obtain the final form of eq. ( |A.12| ) we stih have to perform a cyclic shift by —i places in the indices 
permuted by 7. Using the second form of eq. ( A.16| ) we obtain 



^_^p + ,(,_,) ^ ^_^^M^_ 1(1+1) +^^t;^0-) + ^(.-fc) 

= (_l)E-::o-+7(i)) ^ (A.17) 

which leads us from eq. ( A.15| ) to eq. ( A.12| ). 



B Some lemmas for Consistency Condition II 

Lemma 3: Let A be an n x n matrix with elements ati and Aj be the same matrix with the j'-th row 
shifted by one unit to the left: {Aj)ki = ati for k ^ j and {Aj)ji = aji^i. Expanding the determinant 
of A with respect to the last column we have 

n 

detA = Y^i-ir+Uir^detAl (B.l) 

1=1 

where A^^ is the algebraic complement of a;„. It then holds 

n n 

Y^detA, = J2{-ir+^a,n+idetA*^ (B.2) 

i=i i=i 

Proof: We have to show that when summing over the determinants of the shifted matrices Aj and 
expanding each of them with respect to the last column only determinants of the algebraic complement 
of A remain, all the other terms containing determinants of shifted matrices cancel. 

Starting from the definition we obtain after rearranging terms 

J^det^, =^(-ir+%„+idet^*„ + 5:(-ir+^a,„ ^ det{A,)*^ . (B.3) 

j=l j=l i=l j=l ■ jj^i 

Here we have used the same notation for the shifted matrices and applied the fact that {Aj)'!j^ = A*-. 
We will now show by induction that the double sum in eq. ( p.3| ) vanishes. For n = 2 this holds 
trivially. Inside the double sum we sum over shifted matrices of size n — 1 so we can apply the 
statement eq. ( [B.2D for n — 1. Introducing the (n — 2) x (n — 2) submatrix A^„_^.j„ of A, where now 
the last 2 columns and rows i and j are missing we are left with 



E(-i; 



n+i 

"'in 



i=l 



^(-l)'^-i+^,„det4„_i,, + J2 (-l)"+^a,ndet^*„_i. 



(B.4) 



where the split is due to the missing row when expanding twice. It then follows from the obvious 
symmetry A*^^-^^.,,^ = A'j,^_^.^^ that the two sums inside the bracket cancel. 
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Lemma 4: Let A and Af^ be the same as in Lemma 3 and let A^ be the matrix A with the first L 
entries vanishing in the k-th row: {Aj^)ij = for i = 1, . . . ,L, j = k and aij else-wise. Then it holds 
for any fixed L 



fc=i 



det Aj^ = {n- L) det A . 



(B.5) 



Proof: Let us first state the trivial cases apart from L = 0,n. For L = 1 the determinants just differ 
by the first column and can thus be added up, giving (n — 1) det A after pulling out a common factor. 
For L = n — 1 each A^ can be expanded with respect to the row k which then just gives det A. 

We will now proceed by induction. Due to the above remarks the cases n = 2,3 are trivial. If we 
expand the determinants on the left hand side of eq. (|B.5|) with respect to the last column we obtain 



Edet^i 



fe=i 



E E i-ir^'aindetiAi)*^ + Y.{-ir+\rrdetAl 

k=l i=l;ij^k 



4 = 1 



E(-ir+'a4n(n-l-L)detyl*„ + det^ 



(B.6) 



i=l 



In the first step we have used that (^4, 



L\* 



j4*„ and in the second step we have employed induction 



for n — 1. The sum in the last line then gives (n — 1 — L) det A and adds up to the statement eq. ( [B.5| ). 

Lemma 5: Let ^ be a nonsingular (n — 2) x n matrix and b, c, ^ and r] be n- vectors. Then the following 
property holds for determinants of n x n matrices which differ by the last two rows or columns: 



(B.7) 



A 


A 




A 


A 




A 


A 


b 


c 


— 


b 


c 


= 


b 


e 


e 


r] 




r] 


e 




c 


V 



Proof: Using properties of determinants det(i?C) = det i? det C and det^ = det A-^ the statement 
eq. (|B.7|) is equivalent to 



AA^ Ac At] 
bA^ be brj 
iA^ ^c ir] 



AA^ Ac Ai 
bA^ be bi 
r]A rjc r]^ 



AA^ Ai At] 
bA^ bi br] 
cA^ c^ ei] 



(B. 



We will now use extensively that determinants that differ by one row can be added. We obtain for 
the left hand side 



l.h.s. 



AA^ Ae 


Arj 




AA^ 


Ac Arj 




be 





+ 


bA^ 


br] 


— 


^A^ ie 


iri 




^A^ 


ic irj 




AA^ Ae 


Arj 




AA^ 


Ac Arj 




bA^ 


brj 


+ 


bA^ 


br] 


— 





(.V 




^A^ 


ic 




AA^ A^ 


Ac 




AA^ 


Arj Ae 




bA^ bi 





— 


bA^ 


brj 


+ 


r]A^ 







CA^ 








AA^ Ac Ai 

6c 

r/A rje rj^ 

AA^ Ae A^ 

bA^ bi 

?7^ 

AA^ Ae Arj 

bA^ 6?? 

c^ 



AA^ Ac A^ 
bA^ bi 



rjA^ 

AA^ 
bA^ 
rjA^ 



+ 



rjc rj^ 

Ac Ai 

bi 

rje 

AA^ Ai Ae 

bA^ bi 

cry 



:b.9) 
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In the first line the 1. and 3. term cancel after expanding with respect to the last but one row, the 
same happens in the second line expanding the last row. The last line gives the right hand side of eq. 
([B.8D expanded with respect to the last row if we can show that 



AA^ Ai Ar] 
bA^ hi brj 
cA^ 



AA^ Ai Ac 
hA^ hi 



r]A^ 







AA^ Arj Ac 
hA^ hrj 
iA^ 



Performing similar steps as before this can be shown to be equivalent to 



AA^ Ai Arj 
bA^ 
cA^ 



AA^ Ai Ac 
bA^ 
r]A^ 



AA^ At] Ac 
bA^ 

iA^ 



(B.IO) 



(B.ll) 



We will now apply the property det 



B 
D 



C 

E 



det-Bdet(ii^ — DB C) for matrices where B,E 



are quadratic and B is nonsingular. Having E = 02x2 and B = AA^ nonsingular in the upper left 
corner eq. ( |B.11| ) is equivalent to 



det 



bA^ 
cA^ 



(aA^^ ^ {Ai AifA = det 



det 



bA^ 
rjA^ 

bA^ 
iA^ 



(aA^^ ^ {Ai Ac) 
[aA^Y^ [Ar] Ac) 



and thus to 
det ( {Ai Ari) 



bA^ 
cA^ 



det {Ai Ac) 



bA^ 
r]A^ 



det {At] Ac) 



bA^ 
iA^ 



which can be seen to hold by writing out the 2x2 determinants. 



C The partition function for Nf = 4 adjoint rep. fermions 



(B.12) 



We shall here give some technical details on the derivation of eq. (17. 5| ) in the main text. Our starting 



point is the general relation (7.2) with, for equal masses, the matrix A given by eq. (7.4). The first 
observation is that the matrix A, apart from being antisymmetric, also has a mirror symmetry along 
the line perpendicular to (the conventionally defined) diagonal. In other words: 



A_3 

2' 



Al 3 
2'2 



A 



3 1 

2'2 



A 



i 3 

2'2 



(C.l) 



This means that we only have to evaluate 4 independent matrix elements. We choose these to be 
ai2 = A_3 _i,ai3 = A_3 i,ai4 = A_3 a and 023 = A_i 1. (To avoid the cumbersome indices we 

2'2 2'2 2'2 2'2 

have changed the notation from matrix A to matrix a as indicated). 



Our trick is to first Taylor expand all relevant expressions by means of the Taylor expansion of Bessel 
functions, 



h{x) 



b 00 



E 



:h^r 



2'') ^^kT{b + k + l) ' 



(C.2) 
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and then re-express the result for each matrix element as far as possible in terms of Bessel functions 
again, using the same Taylor expansion. The first steps of this procedure give 



ai2 



^ 



2{n+l) 



(n!)2(2n + 3)(n2 + 3n + 2) 



M 



2n+3 



E 

n=0 

1 oo 

-T 

fi^^n\{n + 2)1(271 + 3) 

1 °° rtJ- f2n+2 



fJ' 



n=0 



n!(n + 2)! 



1 r^f^ 

1 1 /"^M 

-/i(2/x) - — / dtIo{t) , 



after using a simple Bessel function identity. We next evaluate 023: 

00 2n 



5„(n!)2(2n + l) 

= — / dtloit) , 

2/i Jo 



which shows that we have the relation 



ai2 + 023 



f^ 



-hip) ■ 



(C.3) 



(C.4) 



(C.5) 



Next, we evaluate 013 using the same procedure: 

(2n + 2);u2n+i 



ai3 



E 



^J(n + l)!)2(2n + 3) 



1 
1 

1 

M 
1 



00 ^2(n+l) 



A^ 



2(n+l) 



2((n + l)!)2 E((^ + i),)2(2„ + 3) 



^ ((n)! 



.71=0 



1^2 



1 + E 



Ai 



2(n+l) 



^J(n + l)!)2(2n + 3)^ 



^o(2m) - E 



M 



2n 



n=0 



;n!)2(2n + l) 



- [/o(2/i) - 023] , 



which shows that also 013 and 023 are simply related. Using the Taylor expansion of ai4, 

^ (6n + 1)^2" 

"'' " ^o(n!)2(2n + 3)(2n + l) ' 

we find that the following identity holds: 

1 
ai4 - 023 + -ai3 = 2ai2 . 



(C.6) 



(C.7) 



(C.8) 
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Proof: We simply note that the Taylor series for the matrix elements ai4, 023 and ais combine in a 
simple way: 



E 



(6n + l)n 



2n 



M 



2n 



1 (2n + 2)^2n+l 



n=0 L 

00 2n 



(n!)2(2n + 3)(2n + l) (n!)2(2n + l) /x ((n + l)!)2(2n + 3) 
6n + 1 1 2 



n=0 



n2 



n! 

00 ,,2n 



+ 



E 



fj- 



(2n + 3)(2n + l) (2n + 1) (n + l)(2n + 3) 



^Q (n!)2 (2n + 3)(2n + l)(n + 1) 

;^Q {n\)'^{n + l)(2n + 3)(n + 2) 
2ai2 , 



(C.9) 



We are now ready to evaluate the Pfaffian of eq. ([7.2|), using the fact the a Pfafiian of a 4 x 4 
antisymmetric matrix is 



which proves the identity eq. (|C.^). 



Pf(^) 



CI-120'34: — 013024 + 014023 

2 2 

'^12 ~ ^13 + O14O23 , 



(CIO) 



where in the second line we have used our additional symmetry ( C.l ). 

It is handy to evaluate instead /i^Pf(A). We then get successively, using relations ( C.5[ ), ( p.6| ) and 



//3pf(A) 



/i-^[ai2 -0^3 + 014023] 

Hlii2tif + /i^a^g - 2/i2/i (2^)023 - t^Wf^f 

+Hlo{2n)a23 + Ai^0230l3 + Oi4023M^ 

fi[Ii{2nf - Ioi2fif] + ///o(2/x)o23 
+023[M^023 - 2^?Ii{2n) + /i^ai3 + ^^014] 
^i[h{2^if - Ioi2fif] + /u/o(2^)o23 
+023 [/^^023 - 2/i^(/i[ai2 + 023]) + /U^oi3 + f/au] 
/i[Il(2^)2 - /o(2^)2] + ^/o(2Ai)o23 

1 



023 + -013] 



+M 023 [oi4 - 2ai2 - 

fi[Ii{2fif - Io{2fif] + M/o(2Ai)o23 

/i[/l(2/i)2-Io(2^)2 



1 /"^M 
+ -Io(2//) / cit/o(t) 

2 Jo 



(C.ll) 



This is the result quoted in the main text. 
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